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Abstract 

We consider fermion masses and mixings in a renormalizable SUSY 50(10) GUT 
with Yukawa couplings of scalar fields in the representation 10 © 120 © 126. We 
investigate a scenario defined by the following assumptions: i) A single large scale in 
the theory, the GUT scale, ii) Small neutrino masses generated by the type I seesaw 
mechanism with negligible type II contributions, iii) A suitable form of spontaneous 
CP breaking which induces hermitian mass matrices for all fermion mass terms of 
the Dirac type. Our assumptions define an 18-parameter scenario for the fermion 
mass matrices for 18 experimentally known observables. Performing a numerical 
analysis, we find excellent fits to all observables in the case of both the normal and 
inverted neutrino mass spectrum. 
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1 Introduction 



The group SO (10) is a favourite candidate for grand unified theories (GUTs) [T] because 
its 16- dimensional irreducible representation (irrep), the spinor representation, contains 
all chiral fermions included in a Standard Model (SM) family plus an additional neutrino 
SM gauge singlet. Moreover, such theories allow for type I [2 J and type II [3] seesaw 
mechanisms (see also [I]) for the light neutrino masses. In the construction of SO (10) 
theories, there are two options [5], either using low-dimensional scalar irreps but accepting 
non-renormalizable terms in the Lagrangian, or one sticks to renormalizable terms, then 
one has to accept high-dimensional scalar irreps according to P, [7j 

16® 16= (10©126) S ©120 AS , (1) 

where the subscripts "S" and "AS" denote, respectively, the symmetric and antisymmetric 
parts of the tensor product. 

In this paper, we deal with the second option. A special renormalizable model is 
the so-called "minimal SUSY SO (10) GUT" (MSGUT) [8j, which uses, for the Yukawa 
couplings, one scalar in the 10 and one in the 126 irrep in order to account for all fermion 
masses and mixings; it contains, in addition, one 210 and one 126 scalar irrep, in order 
to perform the suitable symmetry breakings. This model has built-in the gauge-coupling 
unification of the minimal SUSY extension of the Standard Model (MSSM). Detailed 
studies of this minimal theory have been performed [9], [10], also with small effects of 
the 120-plet [IT] . Though the MSGUT works very well in the fermion sector, there is a 
tension between the scale of the light neutrino masses and the GUT scale. The reason 
is that the natural order of the neutrino masses in GUTs is u 2 /Mqut ~ 1-5 x 1CT 3 
eV, where we have used v ~ 174 GeV for the electroweak scale and a GUT scale of 
^gut ~2x 10 16 GeV. This neutrino mass scale is too low because J Aml tm ~ 0.05 eV 
where Am^ tm is the atmospheric neutrino mass-squared difference — for reviews on the 
status of neutrino masses and mixing see [13j [H] . Studies of the heavy scalar states [12] , 
together with studies of the fermion mass spectrum, have shown that the MSGUT is too 
constrained [15j, and the tension between the scale of light neutrino masses and the GUT 
scale cannot be overcome: if one has a good fit of the fermion masses, which requires a 
seesaw scale below M G ut, then the gauge coupling unification of the MSSM [TBI HZ] is 
spoiled. 

A natural step for supplying additional degrees of freedom to the MSGUT is to add 
the 120-plet of scalars [18J which appears anyway in Eq. (TjQ) — for early works in this 
direction see The disadvantage is that this step adds a considerable 

number of parameters and reduces the predictability of the theory. Adding the 120 leads 
to a resurgence of the type I seesaw mechanism [10], as a consequence of the collapse of 
the seesaw scale with the GUT scale, because type I seesaw allows to enhance the neutrino 
masses through small Yukawa couplings of the 126 pT8 ] l2l"[ [25] ; without the 120, i.e. in 
the MSGUT, this process leads to the contradictions mentioned above. 

The 120 has electrically neutral components only in its four doublets with respect to 
the SM gauge group. These contribute to the Higgs doublets H u of the MSSM, which 

1 We stress that 10 © 120 alone does not give a good fit in the charged fermion sector [25] , 
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are assumed to be the only light scalar degrees of freedom and the only ones which acquire 
VEVs at the electroweak scale. Thus the MSGUT enlarged by the 120 inherits, from 
the MSGUT, the scalar fields responsible for spontaneous symmetry breaking above the 
electroweak scale. In [24] we took this into account by explicitly making the identification 

WR = M GU T = 2 x 10 16 GeV, (2) 

where wr, which defines the seesaw scale, is the vacuum expectation value (VEV) of 
(10, 1, 3) G 126, with the usual notation for multiplets of the Pati-Salam subgroup [2B] 
S77(4) c x SU(2) L x SU(2) R of 50(10). Furthermore, in that paper we reduced the 
number of parameters in the fermion mass matrices by assuming a horizontal symmetry 
and spontaneous CP violation, i.e. real Yukawa couplings, with CP violation stemming 
from the phases of the VEVs. We showed numerically that this scenario can excellently 
reproduce the known fermion masses and mixings. 

Recently, the role of spontaneous CP violation has been upgraded. A "New MSGUT" 
(NMSGUT) was proposed [27], defined by extending the MSGUT by the 120 and spon- 
taneous CP violation. It was shown that the requirement of spontaneous CP violation 
not only has the virtue of reducing the number of parameters of the theory but it has an 
important impact, via threshold effects, on the unification scale as well; it tends to raise 
the unification scale and with it the masses of all heavy multiplets, thereby suppressing 
baryon decay. 

In the present paper we retain Eq. |2| as a reference point. We do not employ any 
horizontal symmetry but we again motivate real Yukawa coupling matrices by spontaneous 
CP violation. However, we assume that it is of a very specific kind: CP is solely violated 
by imaginary VEVs of the 120; the VEVs of the 10 and 126 are assumed to be real. 
In this way, the mass matrices of the down-quarks, up-quarks, charged leptons and the 
neutrino Dirac-mass matrix are hermitian. This scenario was originally proposed in [21], 
its compatibility with sufficiently slow proton decay shown in [22J . However, in J2TJ [22] it 
was assumed that the type II seesaw mechanism is dominating. Since this is incompatible 
with having only one large scale, we have in the present paper type I dominance and 
neglect possible small contributions of type II, suppressed by v 2 /Mqut- Our scenario 
gives an excellent fit to all known fermion masses, mixings and the CKM phase 5ckm, 
as good as the one in [21], though it is of a rather different type. This shows that the 
fermion data do not fix the enlarged MSGUT in a unique way and there is considerable 
freedom in reducing the number of parameters in this theory. 

The paper is organized as follows. In Section [2] we discuss CP-invariant Yukawa 
couplings and lay out our scenario. The method and results of our numerical analysis are 
discussed in Section [31 In Section @] we present the conclusions. Appendix [A] contains 
a small collection of formulas for the SO(10) spinor representation, which is helpful for 
Section [21 
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2 An 50(10) scenario motivated by spontaneous CP 
violation 



Let us a define a transformation 

i^l(x) — > iCtbUx), , . 

CP: , \ , , N 3 

- e (^)^«(£), 

where ipL is a fermionic 16-plet, C is the charge-conjugation matrix, the 45 gauge fields 
are denoted by WJ* 9 (p < q), e(p) = 1 for p = and —1 for p — 1, 2, 3, x — — x), and 
the 7/ pg are signs. No summation is implied in Eq. ([3]). The 45 (hermitian) generators of 
the gauge group 5*0(10) in the fermionic 16 are given by 

l -a pq = % -T p T q (l<p<g<10). (4) 

For a representation of the operators Y p (p — 1, . . . , 10) of the Clifford algebra and useful 
formulas concerning the spinor irrep 16 see Appendix |A] One can easily check that 
the gauge interaction of the fermionic 16-plet is invariant under the transformation ([3]) 
if ESI ESI 

-VpqVpq = Vpq- (5) 

Since 

T T p =C P T p with e p = (-l) 1+p , (6) 

one finds 

Vpq = te q = (-l) P+q - (7) 

Denoting the generators of the Lie algebra so(10) by M pq , we mention that 

M pq ^ Vpq M pq = SM pq S with S = diag (1,-1, 1,-1, 1,-1, 1,-1, 1,-1) (8) 

is the automorphism associated with root reflection, which is the canonical automorphism 
associated with CP. Such an automorphism exists for all compact Lie groups and is the 
reason why any gauge Lagrangian, whether for fermions or scalars, is CP-invariant (TJ 
[29]. 

Now we transfer the CP transformation to the Yukawa couplings given by the La- 
grangian 

C Y = I (H ab i? aL C- x B TpHpipbL + G ab ipJ L C~ 1 B T p T q T r D pqr il; bL + 



Fab ^ T aL C~ l B T p T q T r T s T t A pqTSt ^ bL ) + H.c. (9) 
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The indices a, b denote the family indices, p, q, r, s, t = 1, . . . , 10 are S'O(IO) indices, and 
C is the charge-conjugation matrix. Summation over family and 5*0(10) indices is implied 
in Eq. 0. The matrix B, which ensures 5O(10) invariance, is defined in Eq. f]A6|) . The 
Yukawa coupling matrices have the properties 

Hab = H ba , G ab = —G ba , F ab = F ba . (10) 
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We define CP transformations (no summations implied) 

H p (x) -> £ p H;(x), (11) 

Dpqr(x) > £,p£,q£,rDpqr(%) > (12) 
A pgrsi (x) — > £p£g£r£s£tA* grst (:r), (13) 

for the scalar fields of the irreps 10, 120 and 126, respectively; the latter two are totally 
antisymmetric tensor fields, A is self-dual in addition. 

Now we require invariance of the Lagrangian under the CP transformation given 
by Eqs. fl3]), (TTTj) . (TT21 and (TT5I) . As an example we take the 10 and obtain 

H ab ^80-%^^ ^ -H ab i)\ L BCV p ^H;r bL = (h^BC-^H^) \ (14) 

The equality sign on the right-hand side of the arrow defines the condition of CP in- 
variance: the CP-transformed Yukawa Lagrangian must be identical with its hermitian 
conjugate. Evaluating Eq. (|14p with the help of Eq. (IA6I) . we find a hermitian Yukawa 
coupling matrix. Performing an analogous computation for the 120 and 126 we arrive at 
the conclusion that the CP transformation requires 

H a b = H^ a , Gab = -G* bal F ab = F b * a . (15) 

Together with Eq. (1101) this means that all Yukawa coupling matrices are real. 

In order to obtain a non-trivial CKM phase 5ckm it is necessary to break CP invariance. 
The scenario we envisage was originally proposed in [21] . In the context discussed here, 
we assume that 

• the VEVs of the 10 and 126 are real, 

• CP is spontaneously broken by the VEVs of the 120, 

• this breaking is maximal, i.e., the VEVs of the 120 are imaginary. 

Thus, the mass matrices of the charged fermions and the neutrino Dirac-mass matrix are 
given, respectively, by 

M d = k d H + iK d G + v d F, (16) 
M u = k u H + in u G + v u F } (17) 
M e = k d H + m t G -3v d F, (18) 
M D = k u H + ik d G - 3v u F, (19) 

with 

k d , u , K dA u,D, v d ,u e R, H = H* = H T ,G = G* = —G T , F = F* = F T . (20) 
The mass matrices f|T6l) - f|T9l are hermitian. The light neutrino mass matrix is given by 
M v = M L - MpM^Mn with M L = w L F, M R = w R F, (21) 
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with scalar triplet VEVs Wl and wr. The mass Lagrangian of the "light" fermions reads 

C M = -d R M d d L - u R M u u l - £ R M e £ L + ^v[Cr x M v v L + H.c. (22) 

We finish this section with some remarks. 5*0(10) models have included the so-called 
D-parity |30j, which is a specific involutory SO (10) transformation which uses the branch- 
ing rule 

16 = (4, 2, 1)©(4,1, 2) (23) 

under the Pati-Salam group [26] and exchanges the two Pati-Salam irreps. One can 
combine CP with D-parity and interpret such a transformation as parity in the usual 
sense [29]. Requiring invariance of the theory under this parity, gives the same restrictions 
on the Yukawa coupling matrices as CP alone, since the theory is invariant under D-parity 
anyway. 

In the CP transformation of the 120, Eq. (Tl2l . we could put a minus sign. Then the 
Yukawa coupling matrix would be hermitian and antisymmetric or, equivalently, antisym- 
metric and imaginaryU In that case, real VEVs of the 120 break CP maximally. 

3 The numerical analysis 

As argued in the introduction, with only one large scale, the GUT scale, in the theory, 
we can neglect the type II seesaw contribution in Eq. (j2T|) FI Then a possible phase of wr 
is irrelevant. With 

H' = k d H, G' = K d G, F' = v d F } (24) 
we rewrite the mass matrices as 



M d 


= H' + iG' + F', 


(25) 


M u 


= r H H' + ir u G' + r F F', 


(26) 


M e 


= H' + ir e G' - 3 F', 


(27) 


M D 


= r H H' + ir D G' -3r F F', 


(28) 


My 


= r R M T D F'- x M D . 


(29) 



Without loss of generality we assume H' to be diagonal. Then all redundant parameters 
are removed and we arrive at 12 real parameters in H', G', F' and six real ratios of VEVs. 
Thus our scenario has 18 independent parameters for 18 observables: nine charged-fermion 
masses, three mixing angles and the CP phase <5ckm in the CKM matrix, the atmospheric 
and solar neutrino mass-squared differences Am^ tm and Am^, and three lepton mixing 
angles. 

Eqs. (l25p — (j29p are amenable to a numerical analysis, which will, in particular, yield 
values for r F and r R . If we fix the triplet VEV wr, e.g. by identifying it with the GUT 
scale — see Eq. ()2]), this analysis will also yield definite values for v d and v u because 

Vd = r R w R , v u = r F r R w R . (30) 

2 This is the choice in [2T]. 

3 A quantitative justification will be given in this section. 
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Quarks 


"I'd 




m s 


19.6 ±5.2 


m b 


1063.6 


u 




m c 


210.3273 +£■«» 


m t 


824 3 3 3 +30267.6 
S^4<M.d _ 1476g 6 


Sl2 


0.2243 ±0.0016 


S23 


0.0351 ±0.0013 


S13 


0.0032 ±0.0005 


#CKM 


60° ± 14° 



Leptons 


m e 


n ^8^ +0.0003 

U - dDS0 _o.0003 




75.6715+0-^8 


m T 


1292.2 


Am| 


(7.9 ±0.3) x 10~ 5 


Am atm 


(2.50 t° f 5 ) x 10- 3 


S 2 
6 12 


0.31 ±0.025 


S 2 
6 23 


0.50 ±0.065 


S 2 
6 13 


< 0.0155 



Table 1: Input data at the GUT scale for M GU t = 2 x 10 16 GeV and tan/? = 10. The 
charged-fermion masses are taken from [51], except for the values of m^, m s and m u ; these 
were obtained by taking their low-energy values from [32J and scaling them to Mqut- As 
for Am 2 tm , we use the value obtained in [J4]. We have copied the remaining input from 
Table I in [17]. Charged-fermion masses are in units of MeV, neutrino mass-squared 
differences in eV 2 . We have used the abbreviations su = sin #12, etc. The angles in the 
left table refer to the CKM matrix, in the right table to the PMNS matrix. 

A reasonable condition on these VEVs is given by [23] 

v \ + vl = (r R w R f (l + rj) < v 2 with v = 174 GeV. (31) 

This inequality certainly holds at the electroweak scale. Assuming that it holds approxi- 
mately at the GUT scale as well, we will subject our fit results to this consistency check. 

To find a numerical solution for the parameters in Eqs. (I2"5l) - (l2"9l . we build as usual [T7J 
1231 E3] a x 2 -function for the 18 observables, 

^(p)=f:[ fi{P l~ Vi )\ (32) 

whose input values yi±Syi are given in Tabled! these values refer to an MSSM parameter 
tan/3 = 10. The letter P symbolizes the set of 18 parameters, i.e. the Yukawa couplings 
and ratios of VEVs. The functions fi(P) express our theoretical predictions, as functions 
of the parameter set P, for the observables, obtained from Eqs. fl25l) - fl29l) . As convention 
for the quark and lepton mixing matrix we use that of the Review of Particle Proper- 
ties [32]. The x 2 -f unc tion is minimized analytically with respect to r R . In this way we 
obtain a x 2 -function of the remaining 17 parameters, which is minimized numerically by 
employing the downhill simplex method |33j . 
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Figure 1: x 2 as a function of the leptonic sin' 



2 #i3 for the normal neutrino mass spectrum. 



In the following we will also consider x 2 functions where a specific quantity is pinned 
down to a given value — for previous uses of such a method see e.g. [TTJ [Ml- If we want to 
pin down a quantity u(P), which is independent of the 18 observables, to a value u, we 
add (uj(P) — u}) 2 /(0.01a;) 2 in Eq. (1321) and minimize the thus obtained xt- If u coincides 
with one of the observables yk, the term above is added but at the same time the term 
(fk(P) — yk) 2 /(8yk) 2 has to be removed from the x 2 °f Eq. (1321 . In that way, we can 
study the sensitivity of our scenario to a variation of a quantity uj. 

A fit in the case of normal neutrino mass ordering: We search for a solution in 
the case of the normal ordering mi < m-i < m$ of the neutrino masses (Am|, = m 2 — m 2 , 
Am 2 tm = m 2 — m 2 ). In that case we find an excellent fit with the following properties: 



In the second relation in this equation we have used wr = Mgut — see Eq. (J2J). The 
corresponding values of the matrix elements of H', G', F' are given by 



X 2 = 0.33, ^v 2 + vl = 26.4 GeV, 



(33) 



mi = 1.81 x 10" 3 eV, m 2 = 0.907 x 10~ 2 eV, m 3 = 5.006 x 10" 2 eV. 



H' 



I 0.198268 






-0.896317 
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R 



Figure 2: \ 2 as a function of R = m min / y Am©, where m min = mi for the normal and m 3 
for the inverted neutrino mass spectrum. 







r F 


r u 


re 




r R 


Normal: 


— 71 .5516 


7.16038 


-2.76118 


6.57185 


5842.373 


1.82618 x lO" 16 


Inverted: 


82.0042 


190.194 


6.60808 


-7.22108 


-43303.97 


2.22386 x 10~ 18 



Table 2: The values of the VEV ratios appearing in the mass matrices ( |25l) -(l29l) obtained 
in the fits with normal and inverted neutrino mass ordering. 



F' 



( 
-2.68402 
-1.67616 

3.23879 
5.69390 



2.68402 


0.580726 

5.69390 
19.1109 



1.67616 
-0.580726 


-10.0743 
-33.5646 



(34) 



-10.0743 -33.5646 -44.3762 



where all numerical values are in units of MeV. The fit values of the VEV ratios are listed 
in Table [2J The \ 2 °f the fit practically comes only from two observables: the pull of 
sin 2 #i3 (leptonic mixing angle) is 0.45 and the pull of is 0.31. 

We can ask the question if our scenario makes some predictions. For the best fit we 
find 5pmns = —2.0°. However, this small value misleading, because pinning <5pmns in X 2 
shows that in the large range —60° ^ <5pmns ~ 60° the fit is still very good, with x 2 ~ 1- 
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The worst x 2 is about 15 and occurs around 5pmns ~ 175°, where for instance m b is not 
well reproduced and the leptonic sin 2 9 13 becomes too large. As for sin 2 # 13 , Fig. [1] shows 
that the preferred value is about 0.01. However, we cannot consider this as a prediction 
since in a wide range around this value the \ 2 is still acceptable. Only at very small 
sin 2 #i3 the fit becomes bad, mainly because of rrid, nib and the atmospheric mixing angle 



#23 . The quantity R = mi/^Amg measures how hierarchical a normal neutrino mass 

spectrum is. The x 2 as a function of R is depicted in Fig. [2j We read off that R ~ 0.2 is 
preferred and the x 2 quickly becomes bad for larger R, mainly owing to m^, m& and the 
leptonic sin 2 #13. Also for very small R the fit worsens, for similar reasons as for large R, 
however, not in a dramatic way. On the other hand, in the MSGUT there is a preferred 
range 0.2 < R < 2 and there is a genuine lower bound on R as well [17] . 



With Eq. (1301) and the upper bound on y v\ + v 2 in Eq. (1331 we see that we are allowed 

to raise Wr to wr ~ 6 x Mgut, without violating the inequality fl3TT) . 

It has to be checked that our numerical solution given by Eq. (jM]) and Table [2] respects 
the perturbative regime in the Yukawa sector. Since the procedure has been explained in 
detail in [24], we confine ourselves to the essentials. The two Higgs doublets of the MSSM, 
Hd and H u , have hypercharges +1/2 and —1/2 and VEVs v cos (3 and v sin (3, respectively. 
The corresponding Yukawa coupling matrices are given by Yd = M&j (v cos (3), etc. It turns 
out that the largest Yukawa couplings are (Y u ) 33 ~ (Yd) 33 , where the largest contribution 
comes from H 33 . Using tan/3 = 10, it is given by thH' zz /{v sin /3) ~ —0.48. This confirms 
that the Yukawa couplings are safely in the perturbative regime. 

Finally, we want to estimate the size of type II seesaw contributions to M. v . The 
corresponding mass matrix is given by wiF'/vd = wlF' / (trMgut) — see Eqs. (El), (ED) 
and fl30l) . The largest element in F' is the 33-element. With the value of this element from 
Eq. (l34"j) and tr from Tabled the type II mass matrix contributes at most (1.2 x 10~ 2 ) x 
wl. Since we expect Wl ~ w 2 /Mgut ~ 10 -3 eV, we find the announced suppression with 
respect to type I seesaw contributions. 

A fit for the inverted neutrino mass spectrum: Searching for a fit by imposing 
the inverted ordering 7713 < mi < m 2 of the neutrino masses (Am^ = m\ — m\ as for the 
normal spectrum, but Am 2 tm = m\ — m 2 ), we find a solution which is even better as in 
the case of normal ordering. It has the following properties: 





X 2 = 0.011, yjv% + v 2 u = 8.46 GeV, 

mi = 4.920 x 10~ 2 eV, m 2 = 5.000 x 10~ 2 eV, m 3 = 2.18 x 10" 4 eV, 



(35) 



with the matrices 



/ 2.39744 




H' 



33.6387 



V 1127.980 



G 




-2.41722 -2.65793 \ 
2.41722 0.0107775 

2.65793 -0.0107775 j 



(36) 
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F' = 



I -1.05201 

-0.0960901 
\ 0.174940 



-0.0960901 
-14.0343 
26.0245 



0.174940 \ 
26.0245 
-52.9848 / 



where all numerical values are in units of MeV, and the VEV ratios are displayed in 
Table [2j For all practical purposes the fit is perfect and there is no need to give any pull 
values. 

Now we come to the predictions of our scenario in the case of the inverted neutrino mass 
spectrum. Concerning CP violation in neutrino oscillations, our best fit gives <5pmns = 
—107.6°. However, this value has no meaning because x 2 as a function of <5pmns is f° r 
all practical purposes. The same is true for the leptonic quantities sin 2 9 2 3 and sin 2 13 in 
the physically relevant ranges. However, there is a definite prediction for the neutrino mass 



spectrum: hierarchy is strongly preferred — see Fig. [2J When the quantity R = m 3 / J Am 2 

becomes large the fit turns bad; however, there is no clear-cut reason, it is mostly the 
down-quark masses and the top-quark mass which are not well reproduced and the fit 
value of the leptonic sin 2 #13 is around its experimental upper bound. 

In the second relation of Eq. fl35l) we have again used our reference value ([2]). Now 
inequality f[3"TT) is respected for wr < 20 x M G ut, i>e. there is more freedom for wr than 
in the normal case. 

As before, large Yukawa couplings in Y u and Yd are induced by if 33. But now, because 
t\d is so big, a slightly larger coupling is (Yd) 13 — irDG' 13 /(v sin (3) ~ i x 0.66, still in the 
perturbative regime. The discussion of the smallness of type II seesaw contributions to 
Ai u proceeds as for the normal spectrum. 

In Eq. ( 13"B1) the elements G' 23 , F[ 2 and F{ 3 are rather small. This might suggest to set 
them zero, which is achieved by the horizontal Z 2 symmetry ipLi — ► —ipLi, D Pq r —> —Dpqr- 
However, this is untenable because it would lead to vanishing 5ckm- The reason is that 
this horizontal can be combined with the CP transformation of Section [2] to a new 
symmetry CP', under which the vacuum state of our scenario is invariant @ Consequently, 
in that case there is no CP violation [541 [55] and one can show — as it must be in such a 
case — that exp(i5 C KM) = ±1 [55] . 

4 Conclusions 

In this paper we have investigated fermion masses and mixings in an SUSY S'O(IO) 
scenariojj originally proposed in [2T], where the Yukawa coupling matrices of the scalars 
in the irreps 10, 120 and 126 are real and CP violation is induced only by imaginary 
VEVs of the 120. This gives a scenario with 18 real parameters in the fermion mass 
sector. Recent results from the MSGUT require a single heavy 5O(10) breaking scale, 
which is then cogent for the MSGUT extended by the 120 as well. 

There are the following differences between Ref. [21] and the present paper: Firstly, our 
M. v is induced by the type I seesaw mechanism and type II is negligible, whereas in [21] 

4 Under CP', the VEVs in<i,u,i,D do not change sign! 

5 In our analysis SUSY enters only via the input parameters whose values we need at the GUT scale. In 
the evolution of these parameters from the electroweak to the GUT scale we assume the renormalization 
group equations of the MSSM. 
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it was assumed that type II dominates. Secondly, we use a purely numerical method, 
employing the minimization of a x 2 function, whereas Ref. [21] uses an approximate 
semianalytical method. 

We have found excellent fits to fermion masses and mixings for both types of neutrino 
mass spectra. We want to emphasize this in particular for the inverted mass spectrum, for 
which the system of fermion mass matrices in the MSGUT — which has no 120 — does not 
allow an acceptable fit [17], though complex Yukawa couplings and VEVs and contribu- 
tions to M. v from both seesaw types are admitted^] The fits presented in this paper have 
the following features. The diagonal Yukawa coupling matrix H of the 10 is strongly 
hierarchical and is responsible, in the charged-fermion mass spectra, for the hierarchy 
between 2nd and 3rd families. The correct size of the neutrino masses is reproduced 
by a cooperation of two effects: rather large contributions thH + ir^G to the neutrino 
Dirac-mass matrix Mr, from the couplings of the 10 and 120, where G is the Yukawa 
coupling matrix of the 120, and a moderately small coupling matrix F of the 126, which 
enters with its inverse in the type I seesaw formula. The contribution of the 120 to the 
charged-fermion masses and to the CKM matrix is rather small, whereas r^G in Mp in- 
troduces large leptonic mixing angles. Similar features were found in the previous sample 
fit of [23], though the assumptions concerning the fermion mass matrices in that paper 
are quite different from those in the present paper, apart from the use of spontaneous CP 
violation in both scenarios. 

Unfortunately, our scenario is not very predictive. However, it does have one clear-cut 
prediction, namely hierarchy for both the normal and inverted neutrino mass spectrum. 



This is quantified by the observable R in Fig. [2] from where we read off m m i n <C J Am| 



Apparently, extending the MSGUT by the 120 leads to an ambiguous situation con- 
cerning fermion mass matrices: quite different assumptions can result in excellent fits. 
Whether these fits are compatible with the NMSGUT [27], where the VEVs are subject 
to certain relations, remains to be checked. One aspect seems to emerge: spontaneous 
CP violation plays in important role in both the fermion mass matrices [24] and the 
spontaneous breaking [27] of SO (10). 

Acknowledgment: W.G. is grateful to L. Lavoura for useful discussions. 



The MSGUT system has 13 absolute values and 8 phases. 
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A The spinor representation of 50(10) 

A possible representation — on the space Tt of the fivefold tensor product of C 2 — for the 
Clifford algebra associated with the Lie algebra so(10) is given by [35] 

IVi = (r { t j) ® ai ® T 2j = 4 5 ~ j) ® a 2 ® l^'- 1 ) (j = 1, . . . , 5), (Al) 

where a superscript (k) denotes the /c-fold tensor product. The 2x2 matrices cr, (i = 
1,2,3) are the Pauli matrices and 1 denotes the 2x2 unit matrix. It is easy to check 
that the T p (p = 1, . . . , 10) fulfill 

{T p ,T q } = 25 pq l^ and T p = Tj = T; 1 . (A2) 
It is well known that the matrices 

i<7 M = ir p r, l<p<g<10, (A3) 
have precisely the same commutation relations as the basis elements 

(M pq ) jk = 5 pj 5 qk - S qj 5 pk , 1 < p < q < N, (A4) 
of so(10). The a pq generate the spinor irrep 16 of so(10) on the 16-dimensional space 

^(i^ + r^n with r u = 4 5) . (A5) 

Note that Tn anticommutes with all T p (p = 1, . . . , 10). 

For the Yukawa couplings one needs the matrix B and its following properties: 

B = iyyW^, B = B ] = B-\ T P B = i P BY p with £ = (-l) 1+p . (A6) 
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